Effect of the equilibrium pair separation on cluster structures 
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A simple pair potential, which equilibrium pair separation can be varied under a fixed inter- 
action range, has been proposed. The new potential can make both face-centered-cubic(fcc) and 
body-centered-cubic(bcc) structure stable by simply changing one parameter. To investigate the 
general effect of the potential shape on cluster structures, the evolution of cluster structures is cal- 
culated for different equilibrium pair separations. The small size clusters(iV < 25), which adopt 
the polytetrahedra, are almost independent on the details of the potential. For the large size 
clusters(25 < N < 150), the potential with large equilibrium pair separation trends to stable deca- 
hedra and close-packed structure, disordered clusters appear for the potential with small equilibrium 
pair separation, while for the middle range of equilibrium pair separation, the icosahedra are domi- 
nated. 



I. INTRODUCTION 

Nanoclusters, with the exotic physical and chemical 
properties, have brought a big chance for future indus- 
try in many areas, it also challenges our current theo- 
ries and models for traditional condensed matters X Until 
now, many phenomena of nano-system still puzzle physi- 
cists, chemists and material scientists. Undoubtedly, any 
precise description of cluster properties needs the cor- 
rect structural models. To date few direct measurements 
of cluster structures are available experimentally^ and 
much of the current theoretical understanding of cluster 
structures has been derived from atomic-scale molecular- 
dynamics(MD) and Monte Carlo(MC) simulations— 

In computer simulations, finding the most sta- 
ble structure corresponds to search the global mini- 
mum of complicated multi-dimensional potential energy 
surfaccs(PES).- The potential energy surfaces of clus- 
ters are usually represented by an appropriate energetic 
model, such as that based on ab-initio method, tight- 
binding(TB) models, or (Semi)empirical interaction po- 
tentials. The model clusters, described by Lennard- 
Jones(LJ) and Morse potentials, have been studied ex- 
tensively (e.g., Refs.— and references therein). Metal 
clusters are also widely studied using both DFT and 
classical many-body potentials(e.g., Refs.— and ref- 
erences therein). The studies were even extended to 
large molecule systems^£^^^^± 4 -> 1 £>± 6 ->± 7 - and multi- 
matter S ystem.i&±&22i2L22i2&2i>2£>2e; ]y[ any exo ti c struc- 
tures, which is probably forbidden in bulk materials, 
have been reported in the previous studies. The excel- 
lent example includes, the cage structures of carbon^ 
the Icosahedra(Iff) 28 ' 29 and decahedra(-Pj7) 30 ' 31 based 
atomic shell structures of metal clusters, and the recent 
found cage structures of gold clusterS ) 32 ' 33 etc. 

Indubitably the cluster structures are determined by 
the details of the interatomic interaction. One of the 
current focuss is to examine the general structural ef- 
fects of the different contributions to the interaction. 
Doye and Wales found that the Friedel type oscillation 
in atomic potentials can strongly modulate the cluster 



structures . 34 ' 35 ' 36 The original and mechanism of disor- 
dered structure in metal clusters were studied in term 
of the different many-body forms i 37 ' 38 Baletto et al have 
shown how the potential forms affect the crossover size 
between different structural motifs (icosahedra, decahe- 
dra, and truncated octahedra).— Michaelian et al have 
made a comparison of the global minimum structure 
for different many-body potentials* 4 ^, Doye and Wales 
investigated the structural consequence issue for a set 
of Sutton-Chen families of potentials.— Gong and his 
coworkers have studied the relativistic effect on the struc- 
ture of gold clusters, which leads the presences of cage- 
like structures ! 32 ' 33 

For studying the general effects of potential shapes on 
cluster structures, Doye and Wales suggested to consider 
a potential which is simple enough that one can compre- 
hend the effects of any changes to its form. This method 
has been used to investigate the effect of the potential 
range and anisotropy on the cluster structures. Braicr 
et al first made the study on six- and seven-atom clus- 
ters bound by the Morse potential for different interac- 
tion ranged 2 - A similar study was preformed on other 
model potentials. Rey and Gallego have shown how the 
structures and melting behavior of hard-core Yukawa 
clusters changing with the range of attractive Yakawa 
taili^ 3 - Based on the generalized LJ, Morse, and Born- 
Mayer potentials, Amano et al recently investigated the 
structure change with the potential shape. 44 Doye and 
Wales made a systematic search of the PES of Morse 
clusters as a function of the interaction range. They 
found that, decreasing the range results in destabilizing 
strained structuresj 45 ' 46 ' 47 ' 48 The general trend have been 
used to explain the growth sequences of other systems.— 

Despite substantial efforts by many researchers, our 
knowledge about the relationship between cluster struc- 
ture and potential shape is still limited on a few factors, 
specially the effect of the interaction range. Although 
the interaction range does play an important role for 
an interatomic potential, other factors also should not 
be neglected. Obviously the further studies along this 
line are needed. Another important factors could be the 
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equilibrium pair separation(dEPs), which reflect the size 
of atoms. Specially, for a fixed interaction range, how 
the equilibrium pair separation(d_eps) affect the cluster 
structure is worthy to answered. Physically, changing 
of cIeps corresponds to changing the atomic size. The 
importance of the atomic size issue has been shown in 
recent studies on Cqq In present paper, by introducing 
a simple model potential, which has the fixed interaction 
range but variable equilibrium pair separation, or equally 
speaking 'the atomic size', we have studied the effect of 
dEPS on the cluster structures. 

The remainder of the paper is organized as follows. 
The following section describes the new developed model 
potential and computational details. The cluster struc- 
tures for a few selected parameters are presented in Sec. 
III. The conclusions drawn from this work are summa- 
rized in Sec. IV. 



nearest neighbors are contributed to energy in fee due 
to the limited interaction range. It needs to point that 
the Morse potential always make the fee more stable over 
bcc. This is the major difference between the present po- 
tential and Morse one. This potential is also similar to 
the Johnson potential for bcc Iron^ which implies that 
this form may represent some major physics of bcc based 
metals. 

To optimize the cluster structure, we first search the 
global minimum among all the known structures for each 
potential and cluster size^ using the steepest-descent 
method. Then we make further exploration for most sta- 
ble structure with the generalized-simulated-annealing 
algorithm, which has been shown in previous work as a 
powerful and efficient procedure^ The most stable struc- 
ture then can be found among these structures. 



II. COMPUTATIONAL DETAILS 

The potential we used was originally obtained from the 
effective pair potentials, of EAM potential for iron ; 51 i 52 
then we re-parameterize it in the current form: 



J £[-£ -2.5] 3 [ 7 - 1.44719(^)]; r« < 2.5^ 



0; Tij > 2.5cr 



where e and a is energy and length unit respectively. 
Tij denotes the distance between atom i and j. 7 is an 
adjustable parameter, which determines the cIeps of the 
potential. The higher value of 7 gives a shorter cIeps- e 
is chosen to keep the potential well depth equal to one 
for each 7, and a equals to one for all case studied. Fig- 
ure 1 shows the new potential with a few selected 7 of 
0.8,0.9,0.95,1.0,1.05,1.1,1,2. For comparing, this figure 
also shows the Morse and LJ potential, where they have 
been fitted to have the same curvature at the bottom of 
the potential well as the present one. From this figure, 
one can see that, changing 7 corresponds to the alterna- 
tion of the equilibrium pair separation, also can be re- 
garded as changing the atomic size. Of course, by chang- 
ing 7, the stiffness of the potential is also changed, which 
is similar to the Morse potential in this point. Compar- 
ing to Morse and LJ potential, the new one is softer in 
repulsive part and stiffer in attractive part, and shorter 
in interaction range. 

Another feature of the current potential is that it can 
make both fee and bcc stable by varying 7. Table 1 
presents the cohesive energy of fee and bcc phase as a 
function of 7. The both fee and bcc have the same en- 
ergy at zero temperature at 7 =1. For 7 larger than one, 
bcc phase is stable over fee, while the smaller 7 favors 
fee phase. The reason is that, for large 7, the poten- 
tial well becomes more and more flat. In this case, the 
dominated interaction in bcc is the both first and second 
nearest neighbors totally 14 atoms, while only 12 first 



III. RESULTS AND DISCUSSIONS 

Fig. 2 plots the second difference of energy(A 2 i? = 
E(N + 1) + E(N - 1) - 2E(N)) as a function of clus- 
ter sizes for all studied 7. From top to bottom, it cor- 
responds to the 7 from 0.8 to 1.2 respectively. For size 
smaller than 24, all the curves have the same trend. This 
implies that the structures are less dependent on the de- 
tails of potentials at small sizes, which is similar to the 
case in Morse potentials^ For all 7 expect 0.8, around 
size from 135 to 147, A2E is almost equal to zero. This 
is because these clusters are based on /f/147 cluster by 
taking a few atoms among its 12 vertex atoms. These 
vertex atoms are all equal due to the symmetry. As the 
size larger than 24, the differences of these curves begin 
to emerge. The curve of 7=0.8 is evidently different from 
others for size larger than 23. 7=0.9 and 0.95 are simi- 
lar each other in whole range of sizes. The major trends 
of 7=1.0, 1.05, 1.1 and 1.2 are close in the most ranges. 
These differences are the reflection of different structure 
type as discussing below. 

Peaks in A 2 E correspond to clusters which are stable 
compared to adjacent sizes and have been found to cor- 
related with magic numbers in mass spectra of clusters. 
We note that there are two types of peaks in A 2 E. One 
type of peaks corresponds to especially stable clusters, 
such as N=13,55,147, the closed-shell IH structures. An- 
other one is directly related to the change of structural 
types. For example, A 2 E for 7=0.8 has several this type 
of peaks for N>100, which actually correspond to the 
close-packed (CP) clusters changing to the DH motifs. 

A few of structural types, i.e, polytetrahe- 
dra(PX'), polytetrahedral involve an ordered array 
of disclinations(PT-rf) IH, DH, CP and disordered(£>/S) 
structures, were found in present study. About these 
structures, the detailed description can be found else- 
where, (e.g., Refs^ and references therein) we briefly 
summary here. PT is made by packing five smallest 
tetrahedra sharing a common edge. The PT can 
further reduce its strain by involve an ordered array of 
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disclinations (hereafter labeled as PT-d), where there are 
more than five smallest tetrahedra sharing a common 
edge. Both PT and PT-d can be naturally divided 
up into tetrahedra with atoms at their corners. IH 
can be decomposed into 20 fee tetrahedra sharing a 
common vertex in the central site. A close-shell IH 
has 20 triangular fcc(lll) faces and 12 vertices. A 
decahedron is made up of two pentagonal pyramids 
sharing a common basis. It has a single fivefold axis 
and is formed by five tetrahedra sharing a common 
edge along the fivefold axis. Both IH and DH can be 
have close geometry shell structure. CP is a pieces of 
fee structure. The truncated octahedron is one of CP. 
By contrast, close-packed structures are composed of 
octahedra and tetrahedra. DIS is that the clusters can 
hard be identified with any order structure. 

Most PT clusters are observed at small sizes. Specially 
all clusters less than 24 atoms are PT type, and most of 
them are identical for all 7 at the same size. This means 
that the cluster structures are less relevant with the de- 
tails of potential. Fig. 3 shows the small size PT clusters 
of N=4-23. These clusters are on a PT growth sequence. 
For example, the 7-atom cluster has a structure in the 
shape of a pentagonal bipyramid, which can be viewed as 
packing of five small tetrahedra by introducing some in- 
herent strains. The clusters larger than seven atoms grow 
by introducing more new tetrahedra on its surface. The 
PT growth sequence keeps and leads to the 13-atom clus- 
ter, which is commonly regarded as an icosahedron, how- 
ever, it is also a PT according to Hoare's description, and 
can be view as the packing of 20 small tetrahedra i 56 i 57 
The 19-atom cluster is packed by small tetrahedra on the 
surface of the 13-atom cluster. The clusters grow further 
by packing around the waist of the 19-atom cluster to 
larger clusters. The PT growth sequence is maintained 
and strain is also accumulated. There are a few clusters 
in the size range smaller than 24, which did not follow 
the growth sequence described above, but still has PT 
structure. The exceptions are also listed in the Fig. 3, 
which is N=8, 21,23 for 7=1.2,1.1, and 0.8 respectively. 
For cluster with N=8 and 7=1.2, which is not the nor- 
mal PT structure, is the fraction of a bec structure. N=8 
cluster can also be viewed as a positive-disclination PT 
cluster(only four tetrahedra sharing a common edge), it 
has more positive strain than normal PT. 

Fig. 4 presents the extension of the PT growth sequence 
based on the small size clusters. 34-atom cluster is com- 
plete form of the waist-packing pattern based on the 19- 
atom cluster. The larger PT clusters(35 < N < 55) 
are based on the 34-atom cluster, which is similar to the 
growth mode form 13-atom PT to 19-atom. 

With increasing of cluster size, the positive strain is 
also accumulated rapidly in PT. To reduce the positive 
strain, PT-d was observed. Fig. 5 presents a few selected 
PT-d clusters. This type of clusters has been discussed 
in several papersj 35 ' 47 i 56 PT-d structure is different from 
PT by introducing six tetrahedra share a common edge 
in some sites. Although disclinations is unfavor in local, 



it does result in the decrease of global strain of clus- 
ters. In fact, this structure is similar to the square- 
triangle Frank-Kasper phases usually used in the quasi- 
crystals»^£i 

The IH clusters appear following the PT-d clusters at 
larger sizes. Some selected IH clusters are listed in the 
Fig. 6. Each cluster is given in the side view and the 
top view. There are many discussions in previous pa- 
pers about IH. 3 The IH clusters can be viewed as pack- 
ing atoms based on a 13-atom pentadechedron but not 
a 13-atom Icosahedron. Usually, comparing with PT, IH 
clusters have larger fcc(lll) surface and shell-like struc- 
ture. These IH clusters are stable due to their extreme 
large fcc(lll) surface and less strains comparing to PT. 
For any IH based cluster, it can always form a closed- 
shell icosahedron by packing enough atoms on it. There 
are also a few IH clusters(see Fig. 7), which are interest- 
ing because they are not the fracture of next closed-shell 
icosaheron, but the fracture of the larger icosahedron. 
For example, IH clusters of N=48,50 are not the frac- 
ture of N=55 icosahedron, but the fracture of N=147 
icosahedron. Also clusters with N=89,90(7 = 0.8) and 
N=62,65(7 = 0.9 - 1.0) are the fracture of N=561 icosa- 
heron. These clusters were discussed by Baletto et al 62 
as natural pathway to the growth of larger IH clusters. 

Fig. 8 shows some DH ,CP and DIS structures. Both 
DH and CP clusters are found for the small value of 
7. For simple pair potentials, the smaller DH and CP 
can only exist for very narrow potential, or potential with 
very large dEPS, such as Ceo clusters. CP definitely have 
the same structure with fee crystal, and DH is also more 
close to fee structure over PT and IH. Disorder clus- 
ters appear for larger 7 around 1.2, where the potential 
is much flat around the bottom. Some disordered clus- 
ters are in fact a serious distorted ordered structures. 
For example, distorted IH (N=85, 103,7=1. 2) are listed. 
There are many disordered structures are hard to be rec- 
ognized based on any ordered structures(N=38,7=1.2). 
There are also some interpenetrated clusters, which ex- 
hibit the combination of two type structures 4£ The in- 
terpenetrated clusters were found during the structure 
motif changing. The remarkable feature of these clusters 
is that it exhibits different structural character by look- 
ing from different direction, for examples, the cluster of 
N=59 and 7 = 1.2 is found between PT and IH. 

Fig. 9 plots the zero temperature structure 'phase dia- 
gram' as a function of both cluster size and 7. It can be 
seen that, the smaller 7, the more DH and CP clusters, 
and the less disorder clusters. For 7=0.8, the structure 
is dominated by the DH and CP. And only one DH and 
one CP clusters are found in the potentials of 7 > 0.9. 
The number of PT clusters increase with increasing 
of 7. For different potentials, this growth sequence 
stops at the different sizes(N=24(7=0.8),N=31(7=0.9), 
N=34(7=0.95),N=37( 7 =1.0,1.05,1.1),N=43(7=1.2)). 
PT-d clusters favor the size within 25 to 107. After the 
PT-d motif first emerges at 7=0.95, its number increases 
with the increasing of 7, and reaches the maximum 
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around 7=1.05. As 7 continues to increase, the distri- 
bution of PT-d motif begins to reduce. At 7=1.2, only 
seven PT-d clusters are found. The disordered clusters 
first emerge at the large size between the two magic 
number N=55 and N=147 at 7=1.0, and its distribution 
increase as the increasing of 7. At the value of 7=1.2, it 
is dominated by the disordered cluster. The Iff-based 
clusters usually appear in larger sizes. IH clusters 
dominate for value of 7=0.9 and 0.95, and its number 
decreases for both 7 larger than 0.95 and smaller 0.9. 

Since clusters have non-neglectful surface effect and 
large deformation in contract with its bulk crystal, the 
competition between the deformation(strain) and surface 
energy plays the key role in determining cluster struc- 
tures. According to this consideration, we can give a 
qualitative explanation for the existing of each structural 
type. Small size clusters have very large surface- volume 
ratio, thus the surface effect is dominated. PT clusters 
are abundant at small sizes due to its the lowest energy 
surface, fcc(lll) surface. The dominant geometric effect 
can not be easily affected unless for enough strong atomic 
interaction. For Morse potential^ only an extreme nar- 
row potential can breaks the PT motif. Since strains in 
PT clusters increase rapidly due to continuously pack- 
ing tetrahedra, PT-d replace the PT with size increasing 
due to the partially release of strains. Further increasing 
sizes, both PT and PT-d are more and more unfavor- 
able, thus the IH clusters appear. IH clusters have the 
similar surface as PT and PT-d, but the inner strain is 
much reduced when taking IH arrangement. This is the 
reason why the IH clusters appear after PT-d ones. 

Among IH, DH and CP structure, IH has the lowest 
surface energy and largest deformation energy, CP is op- 
posite, while DH is in the middle. For small 7(narrow 
potential), the structure with large deformation(strain) 
is unfavorable, CP and DH are dominated. The defor- 
mation(strain) can be neglected for largest 7(more flat 
potential), the PT and disordered clusters are favorable. 
For middle value of 7, the PT-d and IH are appeared. 
Generally DIS should have the largest strain energy, thus 
it can only survive for flat potential, this is the case of 
large 7 in present case. The cluster distribution is in 
agreement with the Doye's qualitative principle that de- 
creasing the range of the pair-potential(the width of po- 



tential well) has the effect of destabilizing strained struc- 
tures. It needs to point that a quantitative analysis is 
necessary for detailed understanding. 

The most stable structure of bulk phase is bec for 7 > 
1 , however only one cluster has the character as a fraction 
of bec, which is the one with N=8 and 7=1.2. On the 
contrary, a large number of IH and PT clusters are found 
for 7 > 1 . The main reasons come from the surface effect 
of small cluster, which result in the clusters adopting the 
most close packed fcc(lll) surface. However, there could 
be the inherent competition between fee and bec, the 
existence of DIS and PT-d clusters could be the results 
for 7 > 1. 

IV. SUMMARY 

In this paper, we have studied how the change of equi- 
librium pair separation (or atomic size) for a fixed in- 
teraction range affects the favored structures of atomic 
clusters. To do so, a simple pair potential, which equi- 
librium pair separation can be varied under a fixed in- 
teraction range by changing one of the potential parame- 
ters has been proposed. This potential also can make 
both face-centered-cubic structure and body-centered- 
cubic stable by changing the parameter. The evolution 
of cluster structures are calculated for several sets of pa- 
rameters. Our results show that, the potential with large 
equilibrium pair separation(larger atomic size) favors DH 
and CP structures, disordered clusters appear for the po- 
tential with small equilibrium pair separation, while for 
the middle range of equilibrium pair separation, the IH 
structures are dominated. Present observation conforms 
the Doye and Walse's qualitative principle that decreas- 
ing the range of the pair-potential has the effect of desta- 
bilizing strained structures. 
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TABLE I: The cohesive energy of both fee and bee for different 
sets of 7 and e, where E bcc and Ef cc are the cohesive energy 
for bee and fee structures respectively, and a6cc and bf cc are 
the lattice constants for bee and fee structures respectively. 
For 7 < 1.0, fee phase is more stable, while 7 > 1.0 bec 
structure is more stable. At 7 = 1.0, fee and bec have the 
same cohesive energy. 
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e 


Ebcc 


Efcc 


Q>bcc 


a/cc 


0.80 


98.496 


-5.326 


-6.000 


2.222 


2.778 


0.90 


29.415 


-5.762 


-6.000 


2.036 


2.557 


0.95 


18.760 


-5.848 


-6.001 


1.932 


2.457 


1.00 


12.786 


-6.036 


-6.036 


1.859 


2.361 


1.05 


9.168 


-6.178 


-6.124 


1.793 


2.264 


1.10 


6.844 


-6.292 


-6.252 


1.730 


2.176 


1.20 


4.172 


-6.604 


-6.547 


1.589 


2.021 
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FIG. 1: The new potential(solid line) with a few selected 7 
of 0.8,0.9,0.95,1.0,1.05,1.1,1,2. For comparing, Morse(dashed 
line) and LJ(dotted-solid line) potential are also shown, which 
has been fitted to have the same curvature at the bottom of 
the potential well as the present one. 
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FIG. 2: The second energy difference (A2E) as a function of 
cluster size for all 7. Peaks in A2E correspond to clusters 
which are stable compared to adjacent sizes. 
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FIG. 3: The small size polytetrahedral clusters of N=4-23. 
Most of them are identical for all 7 at the same size. The 
exceptions are N=8,21,23 for 7=1.2,1-1, and 0.8 respectively. 
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FIG. 4: The extension of the polytetrahedral growth sequence 
based on the small size clusters. 
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FIG. 5: A few selected PT-d clusters. 
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FIG. 6: Some selected IH clusters. Each cluster presented in 
the side view and the top view, with a 6-atom fcc(lll) facet 
marked in the side view. 
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FIG. 7: Some IH clusters are a fraction of larger closed- 
shell icosahedrons. The clusters with N=48 and 50 are the 
fracture of N=147 icosahedron; and N=89, 90(7 = 0.8) and 
N=62,65(7 = 0.9 - 1.0) are the fracture of N=561 icosahe- 
dron. 
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FIG. 8: Some DH, CP and DIS clusters, Top row: CP clus- 
ters; Middle row: DH clusters; Bottom row: DIS clusters. 
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FIG. 9: The zero temperature structural 'phase diagram' as 
a function of both size and 7. From left, 'Filled triangles' for 
PT, 'Stars' for PT-d and 'Vertical lines' for IH. 'Filled square' 
for CP, 'Open diamond' DH, 'Plus' for some minor struc- 
ture including interpenetrated clusters, a truncate tetrahe- 
dra(N=91, 7=0.8)£ and a rare tetrahedron(N=98, 7=0.8).^ 



